Bulk-mediated diffusion on a planar surface: full solution 
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We consider the effective surface motion of a particle that intermittently unbinds from a planar 
surface and performs bulk excursions. Based on a random walk approach we derive the diffusion 
equations for surface and bulk diffusion including the surface-bulk coupling. From these exact 
dynamic equations we analytically obtain the propagator of the effective surface motion. This 
approach allows us to deduce a superdiffusive, Cauchy-type behavior on the surface, together with 
exact cutoffs limiting the Cauchy form. Moreover we study the long-time dynamics for the surface 
motion. 
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I. INTRODUCTION 

Interfaces and the interaction of particles with them 
play a crucial role on small scales in biology and tech- 
nology. For instance, biopolymers such as proteins or 
enzymes diffusing in biological cells intermittently bind 
to cellular membranes, or individual bacteria forming a 
biofilm on a surface use bulk excursions to efficiently relo- 
cate. Similarly the exchange between a liquid phase with 
a solid surface is an important phenomenon in the self- 
assembly of surface layer films and is a ubiquitous pro- 
cess in emulsions. This bulk-mediated surface diffusion, 
schematically shown in Fig. [TJ was previously analyzed 
in terms of scaling arguments and simulations [3-[3| > and 
was unveiled in field cycling NMR experiments in porous 
glasses [5|. Moreover, effects of bulk-surface interchange 
were reported on proton transport across biological mem- 
branes |6|. Recent studies are concerned with effects of 
bulk-surface exchange on reaction rates in interfacial sys- 
tems Q and with surface diffusion of coppper atoms in 
nanowire fabrication 

The remarkable finding of the bulk mediated surface 
diffusion model is that the effective surface motion is 
characterized by a Cauchy propagator 



n a (r,i) 



2n{r 2 + ct 2 f/ 2 



(1) 



where ~ denotes a scaling property ignoring multiplic- 
ative constants, and c is a dimensional factor. The as- 
sociated stochastic transport is of superdiffusive nature 

HUES, 



(2) 




Here we present a strictly analytical approach to this 
process. Our findings corroborate the previous scaling 



Figure 1: Schematic of the bulk mediated surface diffusion: 
The thinner (black) lines show the motion along the planar 
surface with surface diffusivity D s . Following dissociation 
from the surface, thicker (blue) lines depict excursions into 
the bulk volume, in which the diffusion constant is Db. Even- 
tually, the particle rebinds to the surface. 



results for superdiffusion, however, we also derive the 
cutoffs to this behavior: at sufficiently long distances, 
the Cauchy propagator turns over to a Gaussian wing. 
Moreover at longer times the effective surface diffusion 
becomes subdiffusivc, due to the fact that the particle 
spends less and less time on the surface. Normalized to 
the time-dependent surface coverage the effective surface 
diffusion turns over from superdiffusion to normal diffu- 
sion. 



II. DERIVATION FROM A DISCRETE 
RANDOM WALK PROCESS 

We start with a derivation of the coupling between 
surface and bulk in a discrete random walk process along 
the z coordinate perpendicular to the surface (Fig. [1}. 
Let Ni with i = 1,2,... denote the number of particles at 



2 



site i of this one-dimensional lattice with spacing a. The 
number of particles on the surface at lattice site i — are 
termed Ao . The exchange of particles is possible only via 
nearest neighbor jumps. Each jump event between bulk 
sites % — 1,2,... is associated with the typical waiting 
time t. For the exchange between the surface and site 
i = 1 we then have the following law 



at It T dcs 



(3) 



where Td cs is the characteristic time for desorption from 
the surface. The probability of adsorption to the surface 
here is one, i.e., a particle adds to the surface population 
automatically when moving from site 1 to site 0. The 
exchange site at i = 1 between the surface at i = 
and the next bulk site i = 2 is governed by the balance 
relation 

— j^ 1 = A/o - -Ni + —N 2 . 4 

at r dcs t 2t 

Finally, the bulk sites i — 2,3,... are governed by equa- 
tions of the form 

dN 2 (t) 111 
at 2t 2t t 

etc. Let us define the number of "bulk" particles at the 
surface site i = through 



N = —N Q . 

^des 



(6) 



This trick allows us to formulate the exchange equation 
also for site i = 1 in a homogeneous form. Namely, from 
Eq. $5$ we have 



dA/b(t) 1 



(JVi - N ) . 



dt 2t 
Moreover, from Eqs. ([5]) we find 



(7) 



(8) 



for i > 1. 

Let us now take the continuum limit. For that pur- 
pose we make a transition from Ao — > n s as the number 
of surface particles, and JVj —> ant for the bulk concen- 
tration of particles. The factor a may be viewed as the 
lattice constant measuring the distance between success- 
ive sites along the z axis. As the bulk density is a func- 
tion of z while the surface particles are all assembled at 
z = 0, we need this length scale a to match dimension- 
alities between n s and nt, which are [n s ] — cm -2 and 



\n b 



cm 3 . Expansion of the right hand side of Eq. ([7]) 



yields the surface-bulk coupling 



dn s (t) 
dt 



a drib(z,t) 
2r~ a dz 



(9) 



Similarly from Eq. ([8]) we obtain the bulk diffusion equa- 
tion 



dn b (z,t) 
di 



1 2 9 2 Ub 

2> a ~d&' 



Finally, the boundary condition 



—n b (z,t) 



Tde 



(10) 



(11) 



stems from our definition ©. This completes the de- 
scription of the particle exchange between surface and 
bulk, as well as the diffusion of bulk particles along the 
z axis. 

We now turn back to the full three-dimensional prob- 
lem and add to above equations the (x, y) directions (see 
Fig.[T]). As the motion in the perpendicular z direction is 
fully independent, we may simply adjust the Laplacian to 
three dimensions, ending up with the diffusion equation 



dn b (x,y, z,t) 
dt 



dx 2 



d_ 2 

dy 2 



D b -IS + «75 + ^ • ( 12 ) 



If 

dz 2 



with diffusivity D b . For the surface diffusion along the x 
and y coordinates with z = 0, we end up with the two- 
dimensional diffusion equation with diffusivity D s , plus 
the bulk-surface exchange term: 



dn s (x,y,t) 
dt 



s \dx 2 dy 2 



dn b 

dz 



(13) 



z=0 



These two equations are valid in the range < z < oo 
and — oo < x, y < oo, and are supplemented by the initial 
condition 



n a (x,y,t) 



t=o 



n 8(x)6(y) 



(14) 



indicating that initially the particles are all concentrated 
on the surface at x = y = 0. Moreover we observe the 
boundary condition for n b , 



z=0 



n b (x,y,z,t) 
and the bulk initial condition 

n b (x,y,z,t) 



tm s (x,y,t), 



0. 



t=o 



(15) 



(16) 



In what follows, for simplicity of notation we use a unit 
initial concentration uq. 

The surface and bulk diffusivities can be expressed in 
terms of the lattice constant a and the typical waiting 
times between jumps in the bulk, r, and on the surface, 
t s , through 



and 



Dh 



z=0 



a 

6? 



a 



(17) 



(18) 
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respectively. In the continuum limit, both a and r (or 
t s ) tend to zero such that the diffusion constants remain 
finite. In many physical systems bulk diffusion is con- 
siderably faster, such that D s <C Db- Above we also 
introduced the coupling parameter 



DbTdcs 



(19) 



of physical dimension \p] = 1/cm. Small values of /i at 
fixed a and Db correspond to slow bulk-surface exchange. 

For consistency we derive the overall number of 
particles. To this end we define the number particles 
on the surface, N s (t), and in the bulk, Nb(t) through the 
relations 



N s (t) = / dx dyn s (x,y,t) 



and 



poo />oo />oo 

Nf,(t) = / dz dx / dyn s (x,y,t). 

JO J —oo J — oo 



From integration of Eq. (JT3J) we find 

dy n b (x,y,z,t) 



dz 



dt 



— OO </ — OC 



(20) 



(21) 



(22) 



z=0 



Similarly, Eq. JT2J) yields 



dN b (t) 
dt 



D h i dz—^ / dx dyn b (x,y,z,t) 

JO V z J-oo J-oo 

/•OO fOO 



oo J — oo 



Db ~dz ' dx ' d y nb ( x ^y^ z ^) 



(23) 



z=0 



Combination of Eqs. (1221) and (|23|) produces 

j t (N s (t)+N b (t))=0. (24) 



Thus the overall number of particles is conserved, as it 
should be. 



III. SOLUTION OF THE COUPLED DIFFUSION 
PROBLEM 



To solve the set of coupled equations ([11?)) and (fi"2")) 
for the specified initial and boundary value problem, we 
start by defining the two-dimensional Green's function 



G s (x,x',y,y',t) = G s (x -x',y- y',t) 



1 



exp 



( X - X 'Y + (y -yiy 

ADA 



Then, the solution of Eq. (fT5|) becomes 

n s (x,y,t) = / dx' / dy'G s (x,x',y,y',t)n s (x',y',t = 0) 



+ J dt' I dx' I dy'G s (x,x',y,y',t-t') ( D b 



dn b (x',y',z,t') 
dz 



z=0 



With initial condition (fT4")) we thus find 

n s (k x , k y ,s) = G s {k x ,k y ,s) + G s (k x , k y ,s)£? \ Di, 
where 

G s {kx-i ky, s) 



dn b (x,y,z,t) 
dz 



;x -> k x ;y -t k y > ;t -t s 



z=0 



s + D b [kl + fc2] 



(25) 



(26) 



(27) 



(28) 



is the Fourier-Laplace transform of the surface Green's function G s . Here and in the following we denote the Laplace 
and Fourier transform of a function by explicit dependence on the image variables, that is, 



f(s)=J?{f(t);t^s} 



f(t)e~ st dt and g(k) = ^{g(x); x -> k} 



g{x)e dx. 



The bulk particle density according to Eq. (|T2"j) is given by the formal expression [TTJ 

(a; - x') 2 + (y- y') 2 + z 2 



n b (x,y,z,t) 



{AirDi 



x3/2 



4Db(t - 1') 



(29) 



(30) 
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From this expression one can indeed show that, despite the factor z, the coupling equation (|T5|) is fulfilled, see 
Appendix [B] Now, Eq. (j!Z7| requires the derivative of expression (|30[) with respect to z, evaluated at z = 0. To find 
that expression, we first differentiate 



drib \i 



dt' J dx' J dy 



,n s (x',y',t') 



(t-t') 5 / 2 

and then calculate its Fourier-Laplace transform 



exp 



(x-x'Y + (y-y'Y + z 
4D 6 (i - f ) 



9n& 
i9z 



--Ti s {kx-! fcyi s)<^ 



' exp [fcg + fcg] t) 



(31) 



(32) 



z=0 



due to the convolution nature of expression (|31|) . The Laplace transform is evaluated by help of the shift theorem, 
yielding 



f exp(-A[fcg + fcg]t) | = ^/J_ 



We thus finally obtain 



t 3 / 2 



t 3/2 



8->S+D 4 [fc|+fc5] 



9n& 
9z 



-^dI /2 yj ' s + D b (k 2 + k 2 ) n s (k x ,k y ,s). 



(33) 



(34) 



r 



Insertion of this equation into expression (|27[) delivers the Inverse Laplace transformation then yields the exact ex- 
solution for the surface density in Fourier-Laplace space, pression 



i}^x : , s) 



G s {k x , ky, s) 



1 + vDl /2 G s (k x , k y , s) x /s + D b [k% + kl] 

(35) 



N s (t) = e x 2t erfc( X t 1/2 ), 



(40) 



where the Fourier-Laplace transform of the Green's func- where we use the complementary error function 
tion G a was defined in Eq. ([28]) . After some transforma- 
tions we arrive at the exact closed form expression 2 r°° 2 

erfc(z) = -= / e~ (2 d£ = l-crf(z). (41) 

At short times t <C x~ 2 ■, this leads to the initial decay 
N s (t) ~ 1 - ^t 1 ' 2 (42) 



T^s (fax i i *5J — I > 

s + D s [kl + **]+ X <Ja + D b [kl + kl] 

(36) 

with the rescaled coupling parameter 



— r»l/ 2 



-Dfc T dc 



(37) 



of dimension [ X ] = 1/sec 1 / 2 . Relation PJ is the main f the number of surface particles, eventually turning into 
result of this work, and we now consider the consequences \ on g ^j me behavior 



to the surface motion effected by the bulk mediation. 

IV. EFFECTIVE SURFACE BEHAVIOR 

We first determine the number of particles on the sur- 
face, 



N s (t) 



(43) 



oo poo 



oo «/ — oo 



N s (t) -- 

whose Laplace transform is 
N a (s) = n(k x ,k y ,s) 



n s (x 1 y,t)dxdy, 



k x =ky=0 S + X S 



1/2- 



(38) 



(39) 



The asymptotic 1 / y7rt decay stems from the returning 
dynamics to the origin z — of a Brownian motion along 
the z coordinate, i.e., it is proportional to the normaliz- 
ation factor of a one-dimensional Brownian motion. The 
additional prefactor x rescales time with respect to the 
efficiency of the surface-bulk exchange. 

We now turn to the surface dynamics, as quantified by 
the effective mean squared displacement along the sur- 



5 



face. In the Laplace domain, 

(r 2 (*)) s - -V^ >fe n(k x ,k y ,s) 



d 2 1 d 
dk 2 k dk 
4£> s 



n(k x , ky, s) 
2xD b 



(44) 



where k — |k| = ,/k 2 +k 2 . From this expression we 



obtain the limiting behaviors at short and long times. 
Thus, we observe that the short time limit t <C 1/x 2 in 
Laplace domain corresponds to s >• x 2 , such that 



<r 2 (*)) f 



4£> s 

„2 



This translates into the asymptotic time evolution 



4£U 1 



(45) 



(46) 



As long as the ratio Db/D s is sufficiently large, there is 
a superdiffusive component (r 2 (i)^ ~ £ 3 / 2 winning over 
the normal surface diffusion proportional to D s . This 
is exactly the famed bulk mediated superdiffusion ori- 
ginally obtained from scaling arguments by Bychuk and 
O'Shaugnessy [1, 2]. Note that this diffusional enhance- 
ment is accompanied by an almost constant number of 
surface particles, compare Eq. (|42)) . 

Conversely, at long times t ^> 1/x 2 ( or s "C X 2 in 
Laplace domain) we find 



<r 2 (*)>. 



AD S 

9 

X s 



2D h 



\s 



3/2 ' 



corresponding to the temporal behavior 



(47) 



(48) 



Somewhat surprisingly, at sufficiently long times the bulk 
contribution to the effective mean squared displacement 
dominates over the surface contribution for arbitrary ra- 
tio Dt>/D s , giving rise to subdiffusive behavior. This 
subdiffusion occurs due to the ongoing loss of surface 
particles into the bulk, see Eq. (|43|). 

Instead of considering the surface mean squared dis- 
placement (r 2 (i)) s we introduce the normalized effective 
surface mean squared displacement 



<r 2 w>: 



1 



N s (t) 



<r 2 «) s . 



(49) 



This quantity can be interpreted as the surface mean 
squared displacement covered by an individual particle 
that effectively stays on the surface and does not fully 
escape to the bulk. At long times this quantity has the 
limiting form 

^ (r^i))™-^. (50) 
The long time diffusion corrected for the number of es- 
caping particles displays normal diffusion, albeit with the 
bulk diffusivity. 

In the following we neglect contributions from the sur- 
face diffusion proportional to D s , in order not to over- 
burden the presentation. The interesting behavior is due 
to the bulk mediation with weight We quantify the 
motion in terms of fractional order moments, before em- 
barking for the surface propagator. 



Fractional order moments 

We now derive an exact expression for the q-th order 
moments (0 < q < 2) 

(\r\ q (t)) s = J \v\«n s {v,t)d 2 v. (51) 

To this end we utilize the following integral on the plane: 



r 



(52) 



where we used polar coordinates r — |r| and k = |k| with the definition 
corresponding to the two-dimensional vectors r = xe x + 
ye y and k = k x e x + k y e y . Thus we identify 



r q = K(q) j (l - cos(k • r) 



dk 



(53) 



K(q) = 2tt 



1 - J (z) 



dz 
7i+i 



?-(?)K'+!) 



(54) 
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With this trick we can rephrase the gth order moment (|51j) as follows, 

(r q (t)) s = K(q) j J (1 - cos (k • r)) n 8 (r, t)-^dv = K(q) J J n s (r, t)dv - / cos (k ■ r) n s (r, t)dr 



dk 



(55) 



The integral over d 2 r of the surface density is but the 
number of surface particles, such that 



(r"{t)) a =K(q) 



d 2 k r 



k 2 + q 



N s {t)-Re{n s (k,t)} , (56) 



where we replaced the Fourier cosine transform of n s (r, t) 
by the real part of the exponential Fourier transform. 



r 



With the Fourier-Laplace transform (|36|) of the surface 
propagator with D s set to zero, for the Laplace transform 
of the gth order moment we obtain 



(r"{s)) s = K(q) / [n s (k = 0,s)-n s (k )S )] 



dk 



2irK(q) 



s + xVs s + xVs + D b k 2 



dk 



2nK{q) 
2nK(q) 

2irK(q) 



yjs + D b k 2 - y/s dk 

xD b r 00 i 



dk 



s + XV* Jo (Vs + y/s + D b k*) (s + xVs + D b k' 2 ) ki- 1 

9/2 pOD ^ 



X 



s + xVs V s 



dy 



7->9/ 2 



(l+<?)/2 («, _ ^2) 



[h{q)-h{q)}. 



(57) 



where on the way we introduced the substitution y 
\J D b /sk. The two integrals Ii are defined by 



h{q) 



00 y l- q 



1+ s/l + y [ 



-.dy 



(58) 



and 



h{q) 



,1-9 



Vs/x + V 1 + y [ 



- T dy. (59) 



Now we analyze the temporal behavior of the gth order 
moment at short and long times. 



1. Long time behavior 
At long times t ^> 1/x 2 ( or s/x 2 <C 1) we see that 

Jo V 1 + V 1 + v 2 V 1 + y 2 J y 



corresponding to the time evolution 



(62) 



Taking normalization by the number of surface particles 
into account, we find 



9 (s)) no!lu - 



'(*)>* 



N s (t) 



t l/2 



(63) 



That is, at long times the surface diffusion exhibits nor- 
mal scaling behavior. 



2. Short time behavior 

The more interesting case is the short time behavior 
corresponding to the limit t <C 1/x 2 (or s/x 2 ^> 1). Here 
we consider three separate cases: 

(i) The case 1 < q < 2: Since both integrals I\ and I2 
converge, we may simply neglect ^{l)- Then 



and 



)/2 : 



(60) 
(61) 



-(3+«)/2 



(64) 



such that after Laplace inversion we find 

(r q (t)) s ~^ +1 )/ 2 . (65) 
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(ii) The case < q < 1: Now we should take into 
account both integrals, 



h{q)-h{q) 



y 1 q dy 



x Jo (i + v/TTF) + 



(66) 



To estimate the main contribution from this difference, 
we split it into three parts, namely 

h{q)-h{q) ~ 

^( t ...dy+ [ S/X ...dy+ r ...dy\ (67) 
X [Jo Ji JV~slx J 

and evaluate each contribution separately. We find 



I 



(68a) 



and then 

-vVx 



Finally, 



i ^/xr 1 vvV s 

(68b) 



dy fx 



Be) 

Thus the main contribution come from Eqs. (|68b|) and 
5cl). and 



7 1 (g)-7 2 (g)~ S ( 1 -«)/ 2 . 



(69) 



With this estimate the Laplace transform of the qth order 
moment, Eq. ([57]). has the leading order behavior 



This corresponds to 



,(i-9)/2 i 



,(3+?)/2 - s l+q ■ 



(r"(t)) a =s f 



(70) 



(71) 



in the time domain. 

(iii) The case g = 1: This special case requires some 
care. We start with the substitution y = tan 6 in 
Eqs. ([5"5f and (f59"|) . Then the integrals h become 



,7T/2 

/i(1)= / 
Jo 



cos 6(1 + cost 



and 



7T/2 



/ COS 0(1 + ^scos0/x)' 
We can now rewrite Eq. (|57[) in the form 



(72) 



(73) 



(r( S )) s = 2^(1)^^ 



1/2 



s(s- x 2 ) 



tt/2 



X Vo 1 + \/s cos 0/x 



1 , (74) 



where we used [lc 



"/ 2 # 

, | 1 + COS ( 



= 1. 



(75) 



For s/x 2 < 1 the integral in the parenthesis of Eq. ([71)) 
becomes uM 



tt/2 



1 + y/s cos 0/x yi - s/x 2 
Thus, for Eq. (J7JJ| we find 



arctan - 



(76) 



2^i 



27T^(1) 



1 



XD 



1/2 



s(s - x 2 ) 



: arctan < 



X ^/l - s/x 2 ~ V 1 + v^/x 
so that we find the s — > behavior 



- 1 



D 



1/2 



(77) 



(78) 



After Laplace inversion, 



D 



1/2 



and (r(t)Y 6 



(A,*) 1 / 2 (79) 



at long times, t 3> l/x 2 j consistent with Eq. Con- 
versely, for s/x 2 > 1 w e employ 13] 



" /2 # 

l + Vscos0/x 



x In 



^s/ X 2 - 1 + y^/x - 1 
v/s/x 2 -l + l-Vi/x' 



Thus we find 
(r(s)) s = 2nK(l) 



xD h 



1/2 



s(s - x 2 ) 



Vs/x 



In 



V s /x 2 - 1 

We thus obtain the limiting form at large s 



s \ X 



Back-transformed this results in 

(r(t)) s « (r"( S ))r 



tint 



(80) 

- 1 
(81) 

(82) 
(83) 



at short times t <C 1 /x 2 • 

Summarizing our results for gth order moments, at 
short times t <C 1/x 2 the effective surface diffusion exhib- 
its anomalous scaling: the gth order moment of the radius 
scales like t q for < q < 1 and t ( i +1 ^ 2 for 1 < q < 2, 
while the first moment includes a logarithmic contribu- 
tion, tint, consistent with the earlier results in Ref. 0. 
At long times t 3> 1/x 2 the gth order moments scale 
normally with time, proportional to t q l 2 . 
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V. SURFACE PROPAGATOR We perform an inverse Laplace transformation along the 

Bromwich path: 

We now turn to the behavior of the surface propagator. 
If we neglect surface diffusion (i.e., D s = 0) the surface 
propagator from Eq. (1361) becomes 

"«(M)= - \ nh2 - (84) 
I 



n 8 (k,i) 



exp(st) ds 



exp(sx 2 t) ds 



Br s + xvs + D b k 2 2iri J Bl s + Vs + ? 2m 

exp(sx 2 t) c?s 



exp(sx 2 t) ds 
r s + \/s — r2 27TZ 



exp(—K 2 x 2 t) I f exp(sx 2 t) ds f exp(sx 2 £) ds 



V« 2 + 1/4 x , i + f _ v « ! 



1 27TZ 



Br 



1 27TI 



exp(— k x t) 
+ 1/4 



1 



s -> i 



1 



6=^/ K 2 + l/4-l/2,t'=tx 2 lV s + fl 

where on the way we introduced the substitutions s' — s/x 2 and k 2 = D b k 2 /x 2 ■ With 

{—}- ; s -> t \ = -L - ae a2 *erfc f aV~t 

's + a I Vtt^ 



a=l/2+- v /K 2 + l/4,t'=tx 2 



85) 



(86) 



we finally arrive at the Fourier transform of the surface propagator, 

f 



n a (k,t) 



2y/n 2 + 1/4 



I + j - ^ j exp 



exp 



tx erfc 



f 



tx 2 erfc 



(87) 



In compact form, 

n s (k,t) = -^—\ae ax2t evic(axVt] 
a + p I V / 

+f3e~ l3x2t eiic (-pxVt) } , (88) 

with 



(89) 



We now calculate the surface propagator n s (r, t) in the 
limits of short and long times. 



A. Short time behavior 

At short times t <C x~ 2 we distinguish between the 
central part of the propagator and its wings. Starting 



with the central part, we thus focus on distances r <C 
(Dbt) 1 / 2 . Due to n 2 = Dbk 2 /x 2 , this means that we 
consider nxVi ^> 1, and therefore k 2 3> 1. Consequently 
we find a « /3 « k, axVi ~ ftxVt ~ K xVt ^> 1, and 
thus 



erfc ( axv^ 



(90) 



exp(-q 2 x 2 t) 

V^ a xVt 



Then, the first term in the curly brackets of Eq. (b ss ) is 
(note that a 2 — a = k 2 ) 



ae ax2 *erfc (axVtj 



1 



(91) 



The error function in the second term of Eq. (|88p is ap- 
proximately 2. Thus, the second term in the curly brack- 
ets of Eq. (|88p prevails, and we obtain 



,(k,t) ~exp (-Kx 2 i) 



(92) 
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The Fourier inversion is accomplished with the help of 
the integral [l3[ 



xe px J (cx)dx = 



( p 2 + c? y/2- 



(93) 



and we obtain 

n s (r,t) 



o 

e~ xVTFbkt J {kr) 



kdk 
~2^~ 



such that 



n s (r,t) 



xD b ' 1 
2ir (r 2 + X 2 D b t 2 ) 3/2 ' 



(94) 



(95) 



This is exactly the two-dimensional Cauchy distribution 
obtained by Bychuk and O'Shaugnessy from scaling argu- 
ments [l| . We see that indeed radius and time are coupled 
linearly in superdiffusive fashion, r ~ t. However, as we 
will see this Cauchy form holds only for the central part 
of the surface propagator. It assumes steeper wings such 
that at all times any moment of the propagator exists. 
We are quantifying these wings below. 

The number of particles populating this cental Cauchy 
domain of the propagator, i.e., within a circle of radius 
(Ai) 1/2 is 



(D b t) 



1/2 



n s (r,t)dr ~ 1 - x 



Vt, 



(96) 



i.e., at sufficiently short times a major portion of the 
particles are contained in the Cauchy domain. Indeed, 
these particles show the superdiffusive mean squared dis- 
placement 



(D„t) 1/2 



r 2 n s (r,t)dr « xD b t 



3/2 



(97) 



consistent with result (|4"o) for D s = 0, the prefactor dif- 
fering by approximately a factor 3/2. 

The outer part of the propagator at short times is 
found for r (Dbt) 1 / 2 . This region can be divided 
into two subregions, namely, k 2 3> 1 corresponding to 
Dbt <C r 2 Db/x 2 , an d n 2 <C 1 corresponding to 
r 2 ^> D b /x 2 - However, it is easy to check that in both 
of these sub- regions we have axVt "C 1 and PxVi <C 1, 
due to the combination of the conditions of short times 
{xVt "C 1) and large radii (nxVt *C 1). Thus both er- 
ror functions in Eq. (|88l> are expanded at small values of 
their arguments, 



2 2 
erfc(±z) -It —t=z ± ^-p 



z 3 + 



(98) 



with z = axVt < 1 and z = flxVi < lj respectively. 
Our result is 

n s (k, t) ~ 1 - -^ X Vt - -^k 2 D bX t 3/2 (99) 



plus terms of order fc 4 and higher. We conclude that at 
the wings the central Cauchy distribution is truncated by 
Gaussian tails whose dispersion grows like i 3 / 2 . 



B. Long time behavior 

This is the limit t ^> 1/x 2 - As this implies axVt ^> 1 
we may expand the error function as 



cr 



rfc(axVi) 



exp(-q 2 x 2 t) 
y/naxVi 



(100) 



We consider the shape of the propagator in the main part 
excluding the region near the origin, i.e, we have 



(101) 



r 1 » DbVt/x » Db/x 2 



In this domain we may use k <C 1, a ~ 1, and n 2 xVi *C 
1. Since /3 = ^/k 2 + 1/4- 1/2 w k 2 , we have /?xV^ < 1, 
and thus the second error function becomes 



cr 



fc(-/3 X V*) « 1. 



(102) 



The surface propagator thus assumes the shape (recall 
that a 2 — a = k 2 ) 



n a (k,t) 



1 / Qx2t exp(-aVi) 
a V ^/iraxVi 
cxp(~K 2 x 2 t) 

CXp(~K 2 X 2 t) 



+ k 2 exp (-/3x 2 



Equivalently, this means that 

cxp(-fc 2 A,t) 



V^xVt 



(103) 



(104) 



Inverse Fourier transform produces the propagator at 
long times, 



n s (r,t) ~ J e-^n s (k,t)^ 



~ -L= r M kr)e- k2D ^. (105) 

From this we obtain the Gaussian form 

n s (M)~^^ZM. (106 ) 



XVTrt 47rD 6 i 

The number of particles contained in this part of the 
long-time propagator follows from 



n s (r,t)dr ~ — \=. 



(107) 



Comparison to Eq. (l4"3"|) shows that these are virtually all 
particles still on the surface. Thus the Gaussian (|106[) 
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indeed dominates the behavior at long times. This state- 
ment is also corroborated from calculation of the surface 
mean squared displacement from Eq. (|106|) in the domain 



> D b Vi/ 



X- 



r^>D„Vt/ X 

4D b t 



r 2 n s (r, t)dr 



xV^t 



1 



(108) 



Again, up to a small correction, this term coincides with 
our exact result for the long time behavior, Eq. 



Our analysis is based on the exact representation of the 
surface propagator in Fourier-Laplace space from which 
systematic expansions are obtained. For the number of 
surface particles we find an initial stagnation, followed 
at longer times by a decrease proportional to 1/Vt. The 
mean squared displacement (r 2 ) s turns from the initial 
superdiffusive behavior proportional to i 3 / 2 to a subdif- 
fusive form (r 2 ) s ~ t 1 / 2 . The latter corresponds to nor- 
mal diffusion, albeit with the bulk diffusivity, if normal- 
ized to the decay of the surface particles. Our analysis 
is complemented with the derivation of fractional order 
moments for the surface propagation. 



VI. CONCLUSIONS 

Complementing our earlier studies of bulk mediated 
surface diffusion on a cylinder [l(| , we presented a strictly 
analytical treatment of the bulk mediated surface diffu- 
sion problem on a flat surface. While our analysis cor- 
roborates previous scaling results of superdiffusion and 
the Cauchy form of the surface propagator, we find cor- 
rections to this behavior sufficiently far away from the 
origin, and at long times. Thus, the Cauchy domain is 
rectified by a Gaussian decay of the extremities of the 
propagator which we derive explicitly. As a consequence 
finite spatial moments of any order remain finite. At long 
times a Gaussian regime takes over, assembling virtually 
all particles. 
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Appendix A: Alternative derivation of the surface 
density 

In this appendix we derive the Levy flight behavior re- 
covered in the main part of the text based on a comple- 
mentary approach developed in Refs. [T3 - [l6j . We start 
with the coupled diffusion equation along the surface 
(equivalent to Eq. (fT3")) ) as 



dn s 
~~dt 



D.(*- 
\dx 2 



dy 2 



T de 



Tde 



dx' 



dy 1 j dt'W hu i k (x-x',y-y',t-t')n s (x',y',t'), (Al) 
o 



r 



where the last term represents rcbinding events of 
particles that previously desorbed from the surface. The 
kernel Whnik(x-,y,t) in the convolution is the probability 
density to attach at (x, y) at time t after desorption from 
the origin at time zero. Together, the last two terms 
of Eq. (IAip represent jumps along the surface mediated 
by the bulk. Note that we have implicitly assumed that 
all particles start at the surface at t = 0. Otherwise an 
additional term would be needed to represent particles 
binding time to the surface for the first time. 

To obtain Wbuik we need to solve the bulk diffusion 
equation, Eq. (fT2|) . with a single particle released from 
the boundary at the origin at time zero. This can be 
achieved by choosing a boundary condition 



Dt 



dn b 
dz 



z=0 



M r dc 



-n b \ z=Q -S(x)S{y)S(t). (A2) 



This equation corresponds to the current of particles 
into the surface, the second term on the right hand 
side representing the initial release. Supplementing this 



by boundary conditions of vanishing density at infinity, 
n b — > as \x\, \y\,z — > oo, and the initial condition of 
finding no particles in the bulk, n b (x, y, z, t = 0) = 0, we 
obtain Wbuik from the solution in terms of the flux into 
the boundary after the initial release: Wbuik(x , y , t) = 
n b (x,y,z = 0,t)/(/iTdes)- 

To solve Eq. (|12p we again Fourier transform along x 
and y (x —¥ k x and y —¥ k y ) and Laplace transform in 
time (t — > s). Together with the initial condition and the 
boundary conditions at infinity it is easily found that 

n b (k x , k y , z, s) = e- z ^+^n b {k x ,k y , z = 0, s), (A3) 



where k 2 = k 2 x + k 2 . Insertion of Eq. (|A3|) into Eq. (|A2 
leads to a bulk density at the boundary which is 



n b (k x ,k y ,z = 0, s) 



1 



(A^dcs)- 1 + D b ^fk 2 + s/D b 

(A4) 
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Multiplication by the rate constant l/(/XTd cs ) then leads 
to 



W hu \ k {k x ,k y ,s) 



1 + [lT des D b y/k 2 + s/D b 



(A5) 



We are now ready to obtain the surface density n s in 
Fourier-Laplace space. From Eq. (|Alj) with the initial 
condition n s (x,y,t = 0) = no6(x)5(y) we obtain 



n s(kxi ky, s) — 



n 



s + D s k 2 + T dc l[l - W hum (k x ,k y , a)] ' 

(A6) 

in which Eq. (|A5[) can be directly inserted. To obtain the 
expression in Eq. (|36[) of the main text one takes the limit 



of the desorption rate constant t^} s and adsorption rate 
constant l/(/^T<j es ) going to infinity while keeping their 
ratio \i fixed, which leads to 



^dcsl 1 ~ W hul]l {k x , ky, s)} -> fiD^k 2 + s/D b . (A7) 



Together with x = fJ-VDb and no = 1 this makes Eq. (|A6[) 
equivalent to Eq. (|36|) . 

As a byproduct of the approach of this appendix we 
remark that we can now obtain the distribution of bulk 
mediated jump lengths. These are given by A(r = 
yfx 2 + y 2 ) = 2-!rrWhuXk(x,y,s = 0). From Eq. flMJ) we 
find [performing the angular part of the inverse Fourier 
transform similarly to Eq. (|95p ] 



A(r) = 2nr 



kdk Jo(kr) 1 7rr 

2n 1 + ^TdcsD b k [iTdcsDb 2([iTdcsDb) 2 



H (1 



/iTdcs-Df) 



Yn 



[iTdcsD b 



(A8) 



where H„ are the Struve functions and Y v are the Bessel functions of the second kind. If we expand A(r) asymptotically 
at large r we obtain a power law 



A(r) 



flTdesDb 



(A9) 



Thus we see that the Cauchy propagator (|95|) . arising for D s = in the short time regime where almost all particles 
are still bound on the surface, is connected with a Levy flight behavior based on a power law jump length distribution. 



Appendix B: Derivation of the coupling equation (|15[) 

We start by rewriting Eq. (|30|) in the form 



f A ^ f\ t A , f , , n s (x\y',t-T) ( ( x - x ') 2 + (y-y') 2 + z 2 \ 

n »^ Z ^ = J^D^ J J J dy ^T 2 6XP { AD^r J ' (B1) 

We know that the Gaussian is a limiting representation of the 8 function: 

lim 1 exp (- ^ X ~ X '^ \ = S(x - x'). (B2) 



r^O y/4-!TD b T V ^D b T 

Let us now introduce the new variable r' = t/z 2 , such that relation (|B1[) turns into 

( A E. f t/Z * 2 A> f A > f a , n s (x',y',t-T'z 2 ) ( (x - x') 2 (y - y') 2 z 2 
n b {x,y,z,t) = . / z dr dx dy , exp 1 



(4nD b ) 3 / 2 J Q J J z 5 (r') 5/2 \ AD b r' z 2 4D b r'z 2 4D b r' z 2 

(B3) 

Rearranging terms, 

( {x-x'f \ ( (y-y'f \ ( L_ 

, [ t/z2 ^n s (x',y\t-r'z 2 ) CXP { AD b r'z 2 CXP { 4D b r'z 2 CXP \ «V 

n b {x,y,z,t) = / dr dx dy - 



j j j (t') 3/2 \J AttDjjt' z 2 y/^D b T'z 2 V^Db 

(B4) 

We are now in the position to take the limit z — > 0, producing (we omit the prime in r') 



exp 



^mrib(x,y,z,t) = J dr j ' dx' j — S(x - x')S(y - y') 



1 



V^D~b 



(B5) 
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which directly reduces to the expression 



exp 



1 



lim n b (x,y, z,t) = (j,n s (x,y,t)-j== J dr ^_ 3/2 b 7 (B6) 

Finally, with the substitution ( = t~ x / 2 with the Jacobian g?£ = —1/ (2t 3 / 2 ), the remaining integral represents the 
normalization of a Gaussian, 



2 



lim n b (x,y,z,t) = fm s (x,y,t)^f= [ d( exp (--£—)= iin s (x,y,t), (B7) 
and the equivalence with the coupling equation (|15p is established. 
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